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ABSTRACT 

We investigated the influence of environment on cluster galaxies by examining the 
alignment of the brightest cluster galaxy (BCG) position angle with respect to the 
host cluster X-ray position angle. The cluster position angles were measured using 
high spatial resolution X-ray data taken from the Chandra ACIS archive, that signifi- 
cantly improved the determination of the cluster shape compared to the conventional 
method of using optical images. Meanwhile, those of the BCGs were measured using 
homogeneous dataset composed of high spatial resolution optical images taken with 
Suprime-Cam mounted on Subaru 8m telescope. 

We found a strong indication of an alignment between the cluster X-ray emission 
and optical light from BCGs, while we see no clear direct correlation between the 
degree of ellipticity of X-ray and optical BCG morphologies, despite the apparent 
alignment of two elliptical structures. We have also investigated possible dependence 
of the position angle alignment on the X-ray morphology of the clusters, and no 
clear trends are found. The fact that no trends are evident regarding frequency or 
degree of the alignment with respect to X-ray morphology may be consistent with an 
interpretation as a lack of dependence on the dynamical status of clusters. 
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1 INTRODUCTION 



It is well established that the major axes of galax y clusters 
tend t o poi n t tow ar d their nearest n e iRhbo u r (e.g. Binggelil 
19821: iFlinI Il987l: iRhee fc KatgertI Il987l: iPlionisI Il994 



West ct al 1995': ' Plionis fc Basilako'j|2002l : I Chambers et al.l 
2000 , 2002 , Hashimoto et al. 2007a). 

Another alignment effect is that between the orienta- 
tion of the brightest cluster galaxy (BCG), or cD galaxy, 



and that of their parent cluster (e.g. Sastrvl[l968l: Dressier! 
1978'; 'Carter fc M etcalfelll980l : lBinggelilll98d : IStrubleHiggoT 
Plionis et al. . ,200^ Si milar alignment is also reported for 



poor groups of ga l axies (iFuUer et al. 19991). Numerical work 



1995 
1997 



Wes t et al.l 
Wes t .199^; 



Il99ll: Iva n Haarlem fc van de Wevgaert| 
Onuora fc Thomas 2000; Splinter ct a,!.: 
Faltenbacher et al.l 20051 ) show that substructure- 



cluster alignments can occur naturally in hierarchical clus- 
tering models of structure formation such as the cold dark 
matter model. 

Unfortunately, all of these previous galaxy-cluster 
alignment studies are using optically-determined cluster 
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position angles, most of them are based on the Palomar 
Observatory Sky Survey (POSS). Despite the importance of 
these optical investigations, individual galaxies may not be 
the best tracers of the shape of a cluster. Problems can arise 
from foreground/background contamination, as well as the 
fact that galaxies contribute discreteness noise. However, 
it is believed that the X-ray emitting gas within a cluster 
traces its gravitational potential (Sarazin 198^). X-ray mor- 
phology may then be one of the best observable phenomenon 
for determining the cluster shape and orientation. Indeed, 
there are several X-ray stud ies for cluster vs. neighbour- 
clust er alignmen t (e.g. 'Ulmer. McMillan fc K owalskl 

1 19891 : (chambers, Melott fc M iller 2000, 2002, Hashimoto 
et al. 2007a), but there are few galaxy-cluster align- 
ment studies using X-r ay mo rpholog y, except for 
iPorter. Schneider fc Hoesseil ([l993) and iRhee fc Latouj 
( 199ll ). where they investigated the BCG-cluster align- 
ment using low spatial resolution X-ray data, as well as 
traditional cluster shape parameter from apparent galaxy 
distribution, and reported a significant alignment. Unfortu- 
nately, previous X-ray studies are mostly based on Einstein 
data. These data are important, but the exposure depths 
are small and the spatial resolution is rather low compared 
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to recently available X-ray data. The low spatial resolution 
may not significantly affect relatively robust measures such 
as position angle in a direct way, but it will critically hinder 
the accurate removal of contaminating point sources and 
the accurate determination of cluster center, and that may 
significantly affect the estimate of cluster X-ray morphology 
including the position angle. Hence, a new investigation 
using deeper X-ray data with much higher spatial resolution 
is needed. 

Here we report a new investigation of galaxy alignment 
with respect to its parent cluster, using the cluster position 
angle and ellipticity determined by high spatial resolution 
X-ray data taken from the Chandra ACIS archive. Mean- 
while, position angle and ellipticity of BCGs are determined 
from optical images taken with Subaru 8m telescope. This 
paper is organized as follows. In Sec. 2, we describe our main 
sample and X-ray measures, in Sec. 3, details of our optical 
data are described. Sec. 4 summarizes our results. Through- 
out the paper, we use Ho = 70 km s~^ Mpc~^, r2„i=0.3, and 
f2A=0.7, unless otherwise stated. 



2 THE X-RAY SAMPLE, X-RAY DATA 

PREPARATION, AND X-RAY MEASURES 

Here we briefly summarize our main X-ray sample. X-ray 
data preparation, and X-ray measures. More detailed de- 
scriptions can be found in Hashimoto et al. (2007b). 

Almost all clusters are selected from flux-limited X-ray 
surveys, and X-ray data are taken from the Chandra ACIS 
archive. A lower limit of z = 0.05 or 0.1 is placed on the 
redshift to ensure that a cluster is observed with sufflcient 
fleld-of-view with ACIS-I or ACIS-S, respectively. The ma- 
jority of our sample comes from the RO SAT Brightest Clus- 
ter Sample (BCS; lEbeling etahl 119981 ) and the Extended 
ROS AT Brightest Cluster Sample (EBCS; lEbeling et all 
l2000h . When combined with EBCS, the BCS clusters rep- 
resent one of the largest and most complete X-ray selected 
cluster samples, which is currently the most frequently ob- 
served by Chandra. To extend our sample to higher red- 
shifts, additional high-z clusters are selected from various 
deep surveys; 10 of these clusters ar e selected from th e 
ROSAT Deep Cluster Survey fRDCS: iRosati et al.|[l99i ). 
10 from the Einstein Ext e nded Medium Sen sitivity Survey 
(EMSS; iGioia et al.lll990l : iHenrv et al.]|l992l). 14 from th e 
160 Square Degrees ROSAT Survev l|Vikhhnin et al.ll 19981 ). 
2 from the Wide An gle ROSAT Pointed Survey (WARPS; 
iPerlman et~al] |2002|). and 1 fro m the North Ecliptic Pole 
sur vey (NEP: ]Henrv et al.l l2006'). RXJ1054 was discovered 
bv iHasinger et all (1991 )7 RXJ1347 was discovered in the 
ROSAT Ml Sky Survey dschindler et al.lll995l). and 3C29 5 
has been mapped with Einstein Henrv fc HenriksenI (|l986l ). 

The resulting sample contains 120 clusters. At the fi- 
nal stage of our data processing, to employ our full analysis, 
we further applied a selection based on the total counts of 
cluster emission, eliminating clusters with very low signal- 
to-noise ratio. Clusters whose center is too close to the edge 
of the ACISCCD are also removed. The resulting final sam- 
ple contains 101 clusters with redshifts between 0.05 - 1.26 
(median z = 0.226), and bolometric luminosity between 1.0 
X 10"" - 1.2 X 10"*^ erg s"^ (median 8.56 x 10** erg s"^). We 
reprocessed the level=l event file retrieved from the archive. 



The data were filtered to include only the standard event 
grades 0,2,3,4,6 and status 0, then multiple pointings were 
merged, if any. We eliminated time intervals of high back- 
ground count rate by performing a 3 cr clipping of the back- 
ground level. We corrected the images for exposure varia- 
tions across the field of view, detector response and telescope 
vignetting. 

We detected point sources using the CIAO routine 
celldetect with a signal-to-noise threshold for source detec- 
tion of three. We removed point sources, except for those 
at the center of the cluster which was mostly the peak 
of the surface brightness distribution rather than a real 
point source. The images were then smoothed with Gaussian 
(7=5". We decided to use isophotal contours to characterize 
an object region, instead of a conventional circular aperture, 
because we did not want to introduce any bias in the shape 
of an object. To define constant metric scale to all clusters, 
we adjusted an extracting threshold in such a way that the 
square root of the detected object area times a constant was 
0.5 Mpc, i.e. const^area = 0.5 Mpc. We chose const =1.5, 
because the isophotal limit of a detected object was best 
represented by this value. 

The morphology of cluster X-ray emission is then char- 
acterized objectively by the position angle, as well as the el- 
lipticity and the asymmetry. The position angle is defined by 
the orientation of the major axis measured east from north. 
Ellipticity is simply defined by the ratio of semi-major (A) 
and semi-minor axis (B) lengths as: 



Elli = 1 - B/A 



(1) 



where A and B are defined by the maximum and minimum 
spatial rms of the object profile along any direction and 
computed from the centered-second moments by the for- 
mula: 



A' = 




(2) 



(3) 



The asymmetry is measured by first rotating a cluster 
image by 180 degrees around the object center, then sub- 
tracting the rotated image from the original unrotated one. 
The residual signals above zero are summed and then nor- 
malized. Please see Hashimoto et al. 2007b for more detailed 
definitions of morphological measures. 



3 OPTICAL DATA 

To determine the position angle and the ellipticity of BCGs, 
we used optical broad band images taken with Supreme-Cam 
(|Mivazaki et al.lll998l ) on the Subaru telescope. The data 
were retrieved from Subaru- Mitaka-Okayama-K iso Archive 
(SMOKA). Reduction software developed by lYagi et ahl 
( j2002| ) was used for flat-flelding, instrumental distortion cor- 
rection, differential refraction, sky subtraction, and stacking. 
The camera covers a 34' x 27' fleld of view with a pixel 
scale of 0'.'202. The photo metry is calibr ated to Vega system 
using Landolt standards (|Landoltlll99^ . We reflne the orig- 
inal astrometry written as WCS keyword in the distributed 



Alignment of Galaxies and Clusters 3 



archival data using using USN0-A2 catalog with positional 
uncertainties less than ~ 0.2 arcsec. The data were taken 
under various seeing conditions, and we used only images 
with less than ~ 1'.'2 seeing. The optical data retrieved from 
SMOKA contains 30 clusters with redshifts between 0.08 - 
0.9, 

Some clusters have observed through many wavebands, 
and that allowed us to investigate the possible variations of 
our measures caused by waveband shifts. We have decided 
to rely primarily on the R band images for this alignment 
study, because we found that the effect of waveband shift is 
negligible. 

The position angle and ellipticity of BCGs are measured 
exactly the same way as the X-ray cluster emission, namely, 
the position angle is defined by the orientation of the major 
axis measured east from north, and the ellipticity is defined 
by the ratio of semi-major and semi-minor axis. Please see 
Hashimoto et al. 2007b (see also [Hashimoto et al.|[l993 ) for 
further details. As a precaution, we investigate the effect 
of superposed small galaxies sometimes lying on top of the 
extended structure of some BCGs, and we found that these 
superposed small galaxies have little effect on our robust 
measures such as, position angle and ellipticity. 



4 RESULTS 
4.1 Systematics 

One of the haunting, yet unfortunately often lightly treated, 
problem of any study comparing complex morphological 
characteristics of astronomical objects is the possible sys- 
tematics introduced by various data quality, exposure times 
and object redshifts. Depending on the sensitivity of mea- 
sures of characteristics, some susceptible measures may be 
seriously affected by these systematics, producing the mis- 
leading results. 

Unfortunately, investigating the systematics on the 
complex characteristics is not an easy task. To investigate 
the systematic effect of, for example, various exposure times, 
one of the standard approaches is to simulate an image with 
a given exposure time by using an exposure-time-scaled and 
noise-added model image. Unfortunately, we need to approx- 
imate the various characteristics of a model to the compli- 
cated characteristics of a real object, (and those character- 
istics are often what we want to investigate) and this is an 
almost impossible task. 

Meanwhile, if we use the real data, instead of the model, 
we will not have this problem. We can simulate lower signal- 
to-noise data caused by a shorter integration time by scaling 
the real data by the exposure time, and adding Poisson noise 
taking each pixel value as the mean for a Poisson distribu- 
tion. However, this simple rescaling and adding-noise pro- 
cess will produce an image containing an excessive amount 
of Poisson noise for a given exposure time, thus lead us to 
underestimate the data quality. This inaccurate estimate of 
noise is caused by the intrinsic noise already presented in the 
initial real data. The intrinsic noise is difficult to be removed 
without sacrificing the fine spatial details of the object. 

Similarly, to investigate the effect of dimming and 
smaller angular size caused by higher redshifts, in addition 
to the rest waveband shift effect, simple rescaling and rebin- 



ning of the real data will not work, because these manipu- 
lations will again produce the incorrect amount of noise. 

Further difficulty associated with simulation using the 
real data comes from the fact that exposure and redshift 
effects are often coupled, because, in the real observation, 
low redshift objects are usually observed with shorter ex- 
posures than high redshift objects. This coupling further 
poses a serious problem, because simple standard method 
of simulating an observation with 'decreased' exposure time 
will force high redshift data to get degraded, which greatly 
reduces signal-to-noise ratio of already low quality high red- 
shift data. 

To circumvent all of these challenging problems, we de- 
veloped a very useful simulating technique employing a se- 
ries of 'adaptive scalings' accompanied by a noise adding 
process applied to the real images. This technique allows us 
to simulated an image of desired fiducial exposure time and 
redshift with correct signal-to-noise ratio without using a 
tricky artificial model image, thus we can easily investigate 
the effect of various image quality and/or easily change the 
real data to common fiducial exposure and redshift for easy 
comparison. Moreover, the technique can provide us with a 
powerful tool for conducting evolutionary studies, enabling 
us to compare the local objects to the high redshift objects 
without degrading photon-expensive high redshift data of 
low signal-to-noise ratio at all. This method is originally de- 
veloped for the comparison of X-ray image data, but can be 
used for almost all kind of imaging data, including optical 
and NIR images. 

Here we briefly describe the method. Please see 
Hashimoto et al. 2007b for further details. To simulate data 
with integration t=tl, an original unsmoothed image (in- 
cluding the background) taken with original integration time 
to was at first rescaled by a factor Ro/(l-Ro), instead of sim- 
ple Ro, where Ro=tl/tO, tO>tl. That is, an intermediate 
scaled image Ii was created from the original unsmoothed 
image lo by: 

- '-(Ay- 

Poisson noise was then added to this rescaled image by 
taking each pixel value as the mean for a Poisson distribution 
and then randomly selecting a new pixel value from that 
distribution. This image was then rescaled again by a factor 
(1-Ro) to produce an image whose signal is scaled by Ro 
relative to the original image, but its noise is approximately 
scaled by y/Ro, assuming that the intrinsic noise initially 
present in the real data is Poissonian. 

Similarly, to simulate the dimming effect by the red- 
shift, an intermediate scaled image Ii is created from the 
background subtracted image lo by a pixel-to-pixel manip- 
ulation: 

i^(^y) = io{^,yrRl 

where 

i?l = [(l + 20)/(l + 2l)]* (6) 

where zO and zl are the original redshift and the new redshift 
of the object, respectively, and B is the background. 

Finally, to simulate the angular-size change due to the 
redshift difference between zO and zl, the original image will 
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Figure 1. Simulating an image with desired exposure time 
and redshift using tlie real data: Even simulating an image 
with prolonged exposure time is possible with our adaptive scal- 
ing method. Here, optical R band images, taJjen with Subaru 
Suprime-Cam, of the BCG at the center of an example cluster 
(Abell 2219) are shown. Images with original and modified expo- 
sure time and redshift are presented with north is up and east is 
left, (a) Original image: exptime(t)=240s, and redshift(z)=0.228, 
(b) Simulated shorter exposure image with t=10s (e) Simulated 
high-z image with z=0.9, t=240s (d) Simulated prolonged expo- 
sure at high-z with t=1092s, z=0.9. 

be rebinned by a factor R2, then intermediate scaled image 
will be created by rescaling the rebinned image by a factor 
1/(R2-1), before the addition of the Poisson noise. For the 
simulation with 'increased' exposure time, this factor can be 
changed to R3/(Ri-R3) where R3 = t2/t0, t2>t0, where t2 
is the increased exposure time, and tO is the original inte- 
gration time, and (R2 -Ra) > 0. The maximum length of 
integration time we can 'increase' (t2maa;) is naturally lim- 
ited by the original exposure time and how much we increase 
the redshift for the redshift-efFect part, and determined by 
the relationship, 

Rl-Rz = 0, (7) 

which is equivalent to the case when no Poisson noise is 
added after the rebinning. Thus, 

fimax = tOi?2- (8) 

This t2max can be also used as a rough estimate of the ef- 
fective image depth. The t2maa; provides an estimate of the 
image depth much more effectively than the conventional 
simple exposure time because t2max is related to a quan- 
tity that is affected both by exposure time and redshift, and 
thus enabling us to quantitatively compare exposure times 
of observations involving targets at different redshifts (e.g. 
100 ksec at z=0.1 and 100 ksec at z=0.9). 

Although we suspected that our ellipticity and position 
angle were quite robust, as a precaution we investigated the 
possible systematics on these measures introduced by vari- 
ous exposure times and redshifts, using our scaling technique 
described above. 

In Fig. 1, we demonstrate our technique of simulat- 
ing desired exposure time and redshift using the real op- 
tical image of the BCG at the cluster center taken with 
Subaru Suprime-Cam. Original and modified exposure time 
and redshift of an example cluster (Abell 2219) are shown 
with north up and east left, where (a) original image: exp- 
time(t)=240s, and redshift(z)=0.228, (b) simulated shorter 




Figure 2. Similarly with Fig. 1, X-ray images from Chandra 
ACIS, of Abell 2219 are shown with original and modified expo- 
sure time and redshift. North is up and east is left, (a) Orig- 
inal: t=41ks, z=0.228 (b) Simulated shorter exposure: t=10ks 
(z=0.228) (c) Simulated High-z image: z=0.9 (t=41ks) (d) Pro- 
longed exposure at High-z: t=188ks, z=0.9. 
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Figure 3. Ellipticity of BCGs is plotted against ellipticity of 
the X-ray morphology of the host clusters. Interestingly, no clear 
correlation is seen. 

exposure image with t=10s, (c) simulated high-z image with 
z=0.9, t=240s, and (d) simulated prolonged exposure at 
high-z with t=1092s, z=0.9. 

Similarly, in Fig. 2, we use the real X-ray images of 
Abell 2219 from Chandra ACIS, and simulated various ex- 
posures and redshifts, where (a) original image with t=41ks, 
z=0.228, (b) simulated shorter exposure; t=10ks (z=0.228), 
(c) simulated High-z image; z=0.9 (t=41ks), and (d) pro- 
longed exposure at High-z; t=188ks, z=0.9. 

Using this technique, we simulated datasets with vari- 
ous exposure times and redshifts, and measured our cluster 
parameters. We found that our X-ray and optical position 
angles are robust against various exposure times and red- 
shifts. Similarly, we found that other morphological mea- 
sures such as, the ellipticity and asymmetry are quite robust, 
as well. 

4.2 Analyses 

Table 1 shows a summary of our optical cluster sample, 
where APA is an acute relative angle between the position 
angle of X-ray (PA_X) and BCG (PA_BCG), namely, the 
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Table 1. Summary of optical cluster sample 


Clustsr 




PA_X 


PA_BCG 


APA'* 


Elli_X 


Elli_BCG 




(redshift) 


^degree} 




(degree J 






a2034 


0.110 


206.9 


22.24 


4.65 


0.15 


0.50 


a2069 


0.114 


327.7 


331.8 


4.18 


0.46 


0.66 


a750 


0.163 


249.0 


249.6 


0.66 


0.14 


0.31 


rxjl720 


0.164 


355.3 


32.08 


36.7 


0.15 


0.34 


a520 


0.203 


192.4 


228.6 


36.2 


0.26 


0.49 


a963 


0.206 


175.9 


347.0 


8.85 


0.15 


0.42 


a2261 


0.224 


225.8 


14.48 


31.3 


0.14 


0.13 


a2219 


0.228 


309.5 


-79.5 


29.0 


0.38 


0.42 


a2390 


0.233 


298.6 


-57.6 


3.76 


0.30 


0.32 


rxj2129 


0.235 


246.0 


65.28 


0.71 


0.21 


0.54 


a2631 


0.278 


258.2 


84.27 


6.07 


0.29 


0.41 


al758 


0.280 


308.7 


85.53 


43.2 


0.47 


0.21 


a2552 


0.299 


201.7 


216.4 


14.7 


0.18 


0.45 


al722 


0.327 


204.4 


357.7 


26.6 


0.27 


0.14 


ZWC13959 


0.351 


333.4 


346.7 


13.3 


0.21 


0.36 


a370 


0.357 


187.2 


86.96 


79.7 


0.37 


0.11 


rxjl532 


0.361 


227.7 


78.24 


30.5 


0.18 


0.37 


ZWC11953 


0.373 


351.2 


306.3 


44.8 


0.24 


0.51 


ZWC12661 


0.382 


159.5 


206.9 


47.3 


0.12 


0.43 


zwcl0024 


0.390 


5,200 


246.6 


61.4 


0.03 


0.48 


i:xj2228 


0.412 


263,7 


56.27 


27.4 


0.21 


0.54 


i:xjl347 


0.451 


359,2 


343.8 


15.3 


0.20 


0.13 


ms0451 


0.540 


279,9 


344.9 


65.0 


0.26 


0.04 


ciooie 


0.541 


228.7 


244.6 


15.9 


0.19 


0.31 


ms2053 


0.583 


304.8 


336.8 


32,0 


0.25 


0.16 


rxjl350 


0.810 


334.2 


308.9 


25.2 


0.20 


0.47 


i:xjl716 


0.813 


236.2 


255.1 


18.9 


0.17 


0.52 


msl054 


0.830 


266,1 


35.96 


50.1 


0.40 


0.51 


i:xj0152 


0.835 


221.6 


231.0 


9.46 


0.58 


0.29 


wgal226 


0.890 


284.5 


264.6 


19.8 


0.09 


0.15 


a: APA is 


an acute anj 


;le of PA.BCG 


PA_X 
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relative position angle differences greater than 90 degree are 
'folded' and changed to be acute ranging between and 
90 degree. Despite the robust nature of our measures, we 
modify, as a precaution, all of the X-ray and optical obser- 
vations to be equivalent to z=0.9 and t=t2max to eliminate 
any possible small systematics, but otherwise to maximize 
the image quality. 

In Fig. 3, the cUipticity of cluster X-ray morphology 
is plotted against the ellipticity of optical morphology of 
BCGs. Interestingly, in spite of expected alignment of two 
elliptical structures, there are a large scatter and we found 
no strong correlation in the relationship between the two 
ellipticitios. 

Fig. 4 shows the acute relative position angle difference 
between cluster X-ray morphology and BCG morphology 
plotted against the ellipticity of BCGs (top panel) and the 
ellipticity of cluster X-ray morphology (bottom panel). Fig. 
4 shows that the position angle difference tends to be smaller 
than 45 degree, implying that BCGs tend to elongated in 
the same direction of the X-ray distribution of their host 
clusters, particularly for clusters exhibiting relatively high 
ellipticity in their optical BCG morphology and/or in their 
cluster X-ray morphology. Meanwhile, for clusters with very 
low BCG or X-ray ellipticity (ellipticity < 0.1) position angle 
are generally poorly determined, and thus position angle 
difference can bo inaccurate. 

Fig. 5 shows the frequency distribution of the position 
angle difference. There is a strong tendency that we have 
more clusters with an angle difference less than 45 degrees, 
consistent with the observation made in Fig. 4. 

To test this trend more rigorously, wo first employed 
the Kolmogorov-Smirnov (K-S) test. The null hypothesis 
here is that our sample can be drawn from a parent pop- 
ulation of random position angle differences. However, the 
K-S test detects the deviation from the parent population 
(here the population of random position angle differences), 
thus it may loose some sensitivity for testing the cluster 
alignment, where it is likely that position angle difference is 
systematically lower than the random sample. To increase 
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Figure 4. The acute position angle difference plotted versus el- 
lipticity of BCGs (top panel) and ellipticity of cluster X-ray emis- 
sion. The position angle difference is determined by the difference 
between the cluster X-ray position angle and the position angle 
of BCG galaxy. Figures illustrate that clusters with the position 
angle difference less than 45 degree tend to be more abundant, 
particularly for high ellipticity BCGs or clusters, implying that 
cluster X-ray emission and optical light from BCG are aligned. 
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Figure 5. The frequency distribution of the position angle dif- 
ference. There is a tendency that we have more clusters with an 
angle difference less than 45 degrees, consistent with the obser- 
vation made in Fig. 4 iiiiijlyiiig that cluster X-ray emission and 
optical light from BCG are aligned. 
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Table 2. Significance levels of the alignment for various luminous 
galaxies 



Statistics 


BCG 


LG2 


LG3 


LG4 


LG5 


K-S 

Rank Sum 


99.93 
99.99 


81.67 
83.06 


49.93 
57.09 


4.22 
54.51 


8.99 
63.33 



0.6 
0.5 
0.4 

X 

E 0.3 



0.0 
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Figure 6. X-ray morphology versus BCG alignment: Cluster X- 
ray asymmetry is plotted against cluster X-ray ellipticity. Large 
solid circles are clusters showing strong alignment with BCG-to- 
cluster position angle difference less than 20 degree, while solid 
triangles are clusters showing the "modest" alignment, but with 
position angle difference between 20 and 45 degree. Crosses rep- 
resent clusters with no sign of the alignment, and small dots are 
clusters in our main sample without optical Subaru data. No clear 
trends are visible regarding frequency or degree of the alignment 
with respect to X-ray morphology, which can be interpreted as 
lack of dependence on the dynamical status of clusters. 



the sensitivity to an alignment signal, as a second test we 
employed the Wilcoxon-Mann- Whitney rank-sum test. The 
null hypothesis of this test is that the position angle differ- 
ence is not systematically smaller or larger than the ran- 
dom sample. Therefore the test is insensitive to an excess 
of angles around the mean (i.e. 45 deg). When applied to 
our sample, both K-S and rank-sum tests show, not surpris- 
ingly the strong alignment signals, and we find that the null 
hypothesis can be rejected with 99.93% and 99.99% confi- 
dence, respectively, thus confirming that BCGs are signifi- 
cantly aligned to the X-ray emissions of the host clusters. 
We have also investigated the alignment of other luminous 
non-BCG galaxies to the X-ray emissions and we found no 
significant alignment. The results are summarized in table 
2, where LG2 is the second brightest galaxy, LG3 is the 
third brightest galaxy, and LGn is the n-th brightest galaxy 
within a projected distance of 1 Mpc from the X-ray center. 

In Fig. 6, we investigated possible dependence of the 
position angle alignment on the X-ray morphology of the 
clusters. In Fig. 6, the cluster X-ray asymmetry is plotted 
against cluster X-ray ellipticity. Large solid circles are clus- 
ters showing strong alignment between the cluster and BCG 
with position angle difference less than 20 degree, while solid 
triangles are clusters showing the alignment, but with posi- 
tion angle difference between 20 and 45 degree. Crosses rep- 
resent clusters with no sign of the alignment, and small dots 



are clusters in our main X-ray sample without optical Sub- 
aru data. No clear trends are evident regarding frequency or 
degree of the alignment with respect to X-ray morphology. 

We have also attempted to investigate possible de- 
pendence of the alignment on cluster redshifts. We found 
that for clusters less than z=0.35, both K-S and rank-sum 
tests show that the null hypothesis can be rejected with 
99.92% and 99.99% confidence, respectively, while for clus- 
ters greater than or equal to z=0.35, alignment signals are 
somewhat weaker that the null hypothesis can be rejected 
with 93.88% and 83.69% confidence, respectively for K-S 
and rank-sum tests. Similarly, we have investigated the de- 
pendence of the alignment on the cluster X-ray bolometric 
luminosity (Lbol), and we did not find any significant trend: 
for clusters with Lbol greater than or equal to 2 x 10''^ erg/s, 
the null hypothesis can be rejected with 93.45% and 96.42% 
confidence, while for clusters with Lbol smaller than 2 x 10''^ 
erg/s, the null hypothesis can be rejected with 99.78% and 
99.90% confidence, respectively for K-S and rank-sum tests. 



5 SUMMARY 

We investigated the influence of environment on cluster 
galaxies by examining the alignment of the BCG position 
angle with respect to the host cluster X-ray position angle. 
The cluster position angles were measured using high spatial 
resolution X-ray data taken from the Chandra ACIS archive, 
that significantly improved the determination of the cluster 
shape compared to the conventional method of using optical 
images. Meanwhile, those of the BCGs were measured using 
high spatial resolution optical images taken with Suprime- 
Cam mounted on Subaru 8m telescope. 

We found a strong indication of an alignment between 
the cluster X-ray emission and optical hght from BCGs, 
while we see no clear direct correlation between the el- 
lipticity of X-ray morphology and optical BCG morphol- 
ogy despite of the apparent alignment of two elliptical 
structures. In the hierarchical structure formation mod- 
els, the alignment effect could be produced by cluster- 
ing models of struct ure formation such as the col d dark 
matte r model (e.g. ISalyador-Sole fc Solaned Il993l: IWestj 



1994 1: lUsami &: Fufimotc 1 19971 : lOnuora fc Thomad 12000. . 
iFaltcnbac her et al.ll2002l . I2OO5I ). The existence of the align- 
ment effects is also consistent with a cosmic struct ure forma- 
tion m odel such as the hot dark matter model (e.g. lZeldovicbl 
Il970h . where clusters and galaxies form by fragmentation in 
already flattened sheet- and fllament-like structures. 

We have also investigated possible dependence of the 
position angle alignment on the X-ray morphology of the 
clusters, and no clear trends are found. If the X-ray mor- 
phology of clusters reflects dynamical status of clusters (e.g. 
Hashimoto ct al. 20o3)i fhe fa-ct that no trends are evident 
regarding frequency or degree of the alignment with respect 
to X-ray morphology may be consistent with an interpreta- 
tion as a lack of dependence of alignment on the dynami- 
cal status of clusters. Primordial galaxy alignments in clus- 
ters can be damped by various mechanisms such as the ex- 
change of angular momentum in gal axy encounters, violent 
relaxation, a nd secondary infall (e.g. iQuinn fc Binnevl[l99^ : 
ICouttj[l99^ ) over a Hubble time. Thus, in highly relaxed 
clusters, we might naively expect to observe weaker primor- 
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dial galcixy alignment, because there has been sufficient time 
to mix the phases. The fact that we do not see any significant 
alignment trend with respect to the X-ray morphology may 
provide an important constrain on these damping scenarios. 
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